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Abstract: We re-consider the gravitino as dark matter in the framework of the Con- 
strained MSSM. We include several recently suggested improvements on: (i) the thermal 
^ production of gravitinos, (ii) the calculation of the hadronic spectrum from NLSP decay 

and (iii) the BBN calculation including stau bound-state effects. In most cases we find an 
upper bound on the reheating temperature Tr < a few x 10 GeV from over-production of 
^i from bound state effects. We also find an upper limit on the stau lifetime of 3 x 10 4 sec, 
which is nearly an order of magnitude larger than the simple limit 5 x 10 3 sec often used to 
avoid the effect of bound-state catalysis. The bound on Tr is relaxed to < a few x 10 8 GeV 
when we use a more conservative bound on 6 Li/ 7 Li, in which case a new region at small 
stau mass at ~ 100 GeV and much longer lifetimes opens up. Such a low stau mass region 
can be easily tested at the LHC. 
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1. Introduction 

Local supersymmetry (SUSY), or supergravity, models predict the existence of a massive 
spin-3/2 particle, the gravitino G, whose mass in general depends on a supersymmetry 
breaking mechanism. Assuming standard big bang cosmology, it was shown early on that 
cosmological constraints require that the mass of gravitino mg be much less than 1 keV, or 
else heavier than some 10 TeV [1,2]. While a primordial gravitino population can efficiently 
be diluted by inflation [3], subsequently the Universe can be repopulated with gravitinos 
via thermal production (TP) processes involving scatterings of Standard Model (SM) and 
SUSY particles in the hot plasma, with the number density proportional to the reheating 
temperature, Tr. If, assuming i?-parity, the gravitino is the lightest supersymmetric parti- 
cle (LSP) and stable, one can also generate gravitinos via non-thermal production (NTP) 
process of freeze-out and decays of the next-to-lightest superpartner (NLSP). Because of 
the gravitino's exceedingly weak couplings to ordinary matter, the latter process usually 
takes place during or after the period of Big Bang Nucleosynthesis (BBN) and involves 
releasing substantial amounts of electromagneticaly and/or hadronically interacting parti- 
cles (the importance of the latter shown to be important in [4, 5, 6]), which could wreck 
havoc to successful predictions of Big Bang Nucleosynthesis (BBN). In order to avoid this, 
and assuming gravitinos to be a dominant component of dark matter (DM), one imposes 
an upper bound on the reheating temperature of Tr < 10 6 ~ 8 GeV [7, 8, 9, 10, 11, 12] (for 
recent updates see, e.g., [13, 5, 14, 15, 16, 6, 17]). 

Considering NTP processes only, constraints on the parameter space of Minimal Su- 
persymmetric Standard Model (MSSM) were derived in [18, 19, 20], basically eliminating 
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the lightest neutralino as the NLSP (unless mg < 1 GeV [16]). In [14] some of us con- 
sidered a combined impact of both thermal and nonthermal production mechanisms in 
the more predictive framework of the Constrained MSSM (CMSSM) [21] and derived both 
the maximum allowed reheating temperature and the viable regions of parameters of the 
model. An improved analysis using full systematic BBN calculation was next conducted 
in [16]. Both showed, in particular, that gravitino dark matter abundance from NTP pro- 
cesses alone [18, 19, 22, 20] can only agree with observations in regions of SUSY parameters 
where a gluino mass is very heavy, in the multi- TeV range, substantially above an approx- 
imate value derived in [23]. In other words, in order to have gravitino DM consistent with 
superpartner masses at or below the TeV scale, a substantial TP contribution to the total 
abundance is necessary. (In [16] it was also shown that the stau NLSP region consistent 
with the cosmological abundance of DM in gravitinos corresponds to a false (local) vac- 
uum of the CMSSM.) Possible solutions to cosmic lithium problems were also investigated 
in [24]. The papers [14, 16] were next followed by similar detailed analyzes [25, 26]. 

More recently it was pointed out [27] that the existence of unstable (with lifetime longer 
than 10 3 seconds), negatively charged electro- weak scale particles alter the predictions for 
lithium and other light element abundances via the formation of metastable bound states 
with nuclei during BBN [27, 28, 29, 30, 31]. These bound state effects were incorporated 
in the gravitino dark matter with charged particle NLSP [32, 33, 34, 35]. 

In light of these recent developments, we present here an updated analysis of the grav- 
itino DM in the Constrained MSSM. Relative to the previous works [14, 16], we make 
improvements in the following aspects: (i) the thermal production of gravitinos including 
the decay allowed by thermal masses [36] as well as scattering processes; (ii) the calculation 
of hadronic and electromagnetic spectrum from NLSP decay including correct implementa- 
tion of left-right stau mixing; (iii) the updated BBN calculation including the bound-state 
effects of charged massive particles. 

As previously in [14, 16], we will treat as a free parameter and will not address 
the question of an underlying supergravity model and SUSY breaking mechanism. We will 
allow to vary over a wide range of values from 0{ TeV) down to the 0{ GeV) scale (as 
most natural in the CMSSM with gravity-mediated SUSY breaking), for which gravitinos 
(at least those produced in thermal production) would remain mostly cold DM relic, but 
will explore at some level also lighter gravitinos, down to 100 MeV. Less massive gravitinos 
could remain DM in some models of gauge mediated supersymmetry breaking [37, 38, 39, 
40]. 

The paper is organized as follows. In section 2 we describe the improvements made 
in the present analysis compared to the previous ones. In section 3 we summarize the 
experimental and cosmological constraints, show our numerical results and derive an upper 
bound on Tr. We summarize our conclusions in section 4. 
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2. Improvements in the analysis 

In this section we list the improvements made in this paper. 
2.1 Gravitino production 

We will consider gravitinos as DM relics which were produced predominantly via the TP 
and NTP mechanisms stated above. 1 

Since the original computation of TP rates in [9] there have been a number of updates 
and improvements [47, 48, 49, 50, 26, 36]. In particular, a gauge invariant computation 
was performed in [50] (with an extension to the three SM gauge groups done in [26]), 
and a different technique to compute gravitino production from decays allowed by thermal 
masses and the effect of the top Yukawa coupling were applied in [36]. In this update, we 
adopt the result of [36], including the three gauge groups, which leads to about a factor 
of two enhancements compared to the previous calculation [49, 50]. For the renormalized 
gauge couplings and gaugino masses at an energy scale Tr, we used a one- loop evolution 
by the renormalization group equation in the MSSM from GUT scale assuming the gauge 
coupling and gaugino mass unification, 

M i (T)=( 9 -^l) 2 m 1/2 

\9GUT J 

The thermal production will therefore depend on the common gaugino mass mi/ 2 and the 
other parameters of the CMSSM: the common scalar mass too, tan/3 and the trilinear soft 
scalar coupling Aq, as well as on the reheating temperature Tr. In our analysis we use an 
expression for f2~ p /i 2 as computed in ref. [36] 

Q^h 2 - 167 ( ^ ) ( Tr ) ( 7{Tr) ) (2 2) 

G 11 -^HlOOGeVv 1 Vl0 10 GcvJ V^/M 2 J ' ^ 

where the gravitino production rate 7 is the sum of three contributions from decay, sub- 
tracted scattering and top Yukawa coupling, 

7 = 7D + 7s" b + 7io P , (2.3) 

and the details are given in ref. [36]. 

Regarding the non-thermal production of gravitinos, we proceed in the usual way. 
Since all the NLSP particles decay after freeze-out, the gravitino relic abundance from 
NTP, f2~ TP /i 2 , is related to ^nlsp^ 2 _ the relic abundance that the NLSP would have had 
if it had remained stable - via a simple mass ratio 

NT P/l 2 = ^^ NLSp/ ^ {2A) 
G "1-NLSP 

1 In addition, there could be other possible ways of populating the Universe with stable gravitinos, e.g. 
via inflaton decay or during preheating [41, 42, 43, 44], or from decays of moduli fields [45]. In some of 
these cases the gravitino production is independent of the reheating temperature and its abundance may 
give the measured dark matter abundance with no ensuing limit on Tr. In general, such processes are, 
however, much more model dependent and not necessarily efficient [46], and will not be considered here. 



(2.1) 
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Note that Q^ TP h 2 grows with the mass of the gravitino vtip. When we calculate the relic 
number density of staus, we treat correctly the mixing of left and right stau states using a 
supersymmetric particle spectrum calculator SuSpect [51] and employ micrOMEGAs [52] 
to calculate a relic density of NLSPs, which might affect the annihilation cross section 
when the maximal mixing happens at large tan/3, as noticed by [35]. 

The total abundance of the LSPs is the sum of both thermal and non-thermal produc- 
tion contributions 

n^h 2 = n% p h 2 + n% TP h 2 . (2.5) 

Ct Ct 

Since it is natural to expect that the LSP makes up most of DM in the Universe, we can 
re-write the above as [53] 

ft d h 2 = njfh 2 (T R , m d ,m- g , m NL SP, • • •) + — ^nlsp^ 2 = ^dm ^ 0.1. (2.6) 

mNLSP 

2.2 Hadronic and electromagnetic spectrum from NLSP decay 

In order to calculate light element abundances during BBN, including the impact of massive 
particle decays, in our case the NLSP, we need to know the detailed decay products and 
their spectrum in addition to NLSP lifetime and relic number density. 

The lifetime of the stau or neutralino NLSP decaying to the gravitino is dominated by 
two-body decays, 

r&t ~r(fi->rG), _ _ (2.7) 

r* ot ~ r( x - g 7 ) + r( x - gz) + r( x -> Gh), (2.8) 

where, for the stau NLSP, 

4 

and for the neutralino NLSP the dominant decay to photon and gravitino is given by [20] 



r(fi -> tG) = -T^-^pr-^T \ 1 — T ) ' ( 2 - 9 ) 
487r M 2 m 2 ~ \ m 2 ~ 1 




^ g = \Nn«*tw + Nudntwf "4 L _ L + . (2 . 10) 

48tt M 2 x m 2 ~ y ml) \ m x) 

These decay widths depend mostly on the LSP and NLSP masses, but also on Nij , the 
neutralino mass mixing matrix. Note that Mp\ denotes the reduced Planck mass. The 
lifetime can be sufficiently long r ~ 10 2 — 10 5 s to be interesting for BBN and the lithium 
anomalies. 

The two-body decays also dominate the electromagnetic cascade with the decay of 
the tau-lepton (producing either directly leptons or mesons that decay electromagnetically 
before interacting strongly with the light elements). The electromagnetic branching ratio 
and energy injected in the cascade for the stau are given by: 

f _ r(f 1 ^rG) f „l( <-rr% \ 

T tot - 1 ' ^"2 I 2m f I [2Ai) 
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The approximate fraction 1 /2 in the total tau energy is a reflection of the fact that some 
energy is lost into non-interacting neutrinos. Similar results are obtained for the neutralino: 

(m\- m 2 ~\ 

B} m - 1, = J ■ (2.12) 

Hadronic energy contributions come from 3-body, or from 4-body, NLSP decays: 

fi -» tGZ/j -» rGgg, (2.13) 
fi -» u T GW -> i/ T Ggg, (2.14) 
fi -» rGTi -> rGgg. (2.15) 

In our previous analyzes [14, 16] we approximated these processes by a simple formula, and 
the process (2.14) was set to zero, since we assumed the stau to be a pure right-handed 
state. Here we incorporate the mixing and use the programs, CalcHEP [54] (an automatic 
matrix element generator) and PYTHIA [55] (a Monte-Carlo high-energy-physics event 
generator) to produce the spectrum of the decay products properly. With this update 
eq. (2.14) is not zero any more due to the left-handed component enabling this process. 

Following [25] we use the qq invariant mass m qq - to calculate the decay width in 
eq. (2.13). For completeness virtual Higgs exchanges were also calculated. The decay 
width is given with a cut on the qq invariant mass m q ^ t = 2 GeV below which no nucleon 
can be produced, 

F(f - rGqq- m%) = d gXf ~> rGqq) 

J m cut dm q q 

The branching ratio is given by 



r(f -> rGqq; m c ^) 

B h (f -> rGqq; m c q f) = -± qq ' ■ (2.17) 

tot 



The hadronic decay is dominated by an exchange of Z/7. The W processes are suppressed 
by one or two orders of magnitude compared to the Z/7 because the left-handed component 
is comparatively small in the CMSSM. The Higgs contribution is even more negligible due 
to the Higgs mass in the propagator. This last contribution can be enhanced at large tan (3 
which increases the Higgs boson couplings but it remains below the W contribution. For 
the neutralino NLSP, the hadronic cascades come from 3-body decay x ~^ Gqq: 

X - GZh - Gqq, (2.18) 
x _> Gh -> Gqq, (2.19) 

and the decay width is calculated similarly as for the stau 

j m cut am q q 
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2.3 BBN calculation 



As in our previous work [16], for each point in supersymmetric parameter space we perform 
a complete BBN calculation of light element abundances with NLSP decay-induced cas- 
cades. These calculations include all thermal nuclear reactions as well as all (nonthermal) 
interactions important for the developments of the electromagnetic and hadronic cascades. 
Details about the effects of the around 100 nonthermal interactions required to attain an 
estimated BBN yield accuracy of 30% may be found in [6] . The calculations are performed 
at n b h 2 = 0.02273, as inferred by WMAP [56]. 

The improvements with respect to ref. [16] are twofold. Firstly, the current BBN 
calculations take full and accurate account of catalytic effects [27, 32, 28, 29, 30, 57, 58, 59, 
31, 60, 61]. Secondly, the energies of injected nucleons are treated in a much improved way. 
It has been recently shown that the formation of bound states between electrically charged 
NLSPs (here principally the stau f) and nuclei towards the end of BBN lead to drastic 
changes of some thermal nuclear rates [27, 30, 62, 61]. Various modifications have been 
proposed, but the by far most important one is the replacement of the weak quadrupole 
transition D + 4 He — > 6 Li + 7 with the catalytic process D + ( He — f) — > 6 Li + r, with a rate 
~ 10 7 times larger than the former. Since both ( 4 He — f ) bound states as well as 6 Li first 
survive photo-disintegration/nuclear destruction at 10 sec, 6 Li overproduction rules out 
most of the stau-NLSP parameter space with stau decay time r^5x 10 3 sec [27, 32, 33]. 
The exception here are very small freeze-out NLSP stau-to-entropy ratios Y = n/s ^ 
10~ 15 — 10 -14 [31] as recently proposed to occur at large tan/? [63]. Other modifications 
to BBN due to the bound state induced catalytic effects are the possibility of appreciable 
9 Be synthesis [60] which has been used to derive limits on abundances of stau NLSPs [62]. 
Though we fully calculate 9 Be/H ratios, we refrain from using them for limits, as the 
catalytic rates seem currently very uncertain [61]. Finally, late-time destruction of 6 Li 
and 7 Li through reactions on (p — f) bound states, as proposed in ref. [59], has proven 
unimportant when the newly determined catalytic rates of ref. [61] are employed. 

In ref. [16] a partition of energy into hadronic three- and four- body decays, such as 
f — > GtZ and f — > Grqq was approximated by Ez ~ E qq ~ (mf — mg)/3. In ref. [25] 
both processes were treated in detail and it was found that the invariant mass-squared 
m qq = {Pq + Pq) 2 °f produced qq flux tubes is peaked around the Z-mass squared m\. It 
was argued that the average of (m q q) should be used as an "effective hadronic energy" for 
BBN calculations rather than E qq « (mf — mg)/3. In particular for heavy mf ~ 1 TeV 
this would make a factor ~ 3 difference in hadronic effects. As described in more detail in 
ref. [64], the situation is actually more complicated. In fact (mqq) is not a good estimate 
of the energy of the qq flux-tube in the cosmic rest-frame. Nevertheless, for independent 
reasons, at early times, r & 300 sec taking (m qq ) is indeed a better approximation, which 
has been adopted in the present analysis. A completely accurate determination would 
require knowledge of the detailed primary nucleon energy spectrum due to the NLSP 
decays. This, is however, numerically currently not feasible for the large number of models 
to be calculated. Some detailed calculations have been performed in ref. [64] indicating a 
~ 30% uncertainty in the results. 
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3. Constraints and results 



Unless otherwise stated, we will follow the analysis and notation of our previous papers [14, 
16] to which we refer the reader for more details. Here we only summarize the main results 
and differences. 

Mass spectra of the CMSSM are determined in terms of the usual four free param- 
eters mentioned above: tan/3, m^, mo, and Aq, as well as on sgn(yu) - the sign of the 
supersymmetric Higgs/higgsino mass parameter [i. The parameter fi 2 is derived from the 
condition of the electroweak symmetry breaking and we take fi > 0. 

We calculate the thermal contribution fi~ p /i 2 following the result of [36]. In evalu- 
ating the non-thermal part 0~ TP /i 2 , we first compute the number density of the NLSP 
after freeze-out (the neutralino or the stau) with high accuracy by numerically solving 
the Boltzmann equation including all (dominant and subdominant) NLSP pair annihila- 
tion and coannihilation channels. For a given value of mg, we then compute 0~ TP /i 2 via 
eq. (2.4). 

After freeze-out from the thermal plasma at t ~ 1CT 12 sec, the NLSPs decay into 
gravitinos at late times which strongly depend on the NLSP composition and mass, on mg 
and on the final states of the NLSP decay [22, 20]. The exact value of the NLSP lifetime in 
the CMSSM further depends on a possible relation between mg and m\/2 and/or mo but 
in the parameter space allowed by other constraints it can vary from > 10 8 sec at smaller 
m-NLSP down to 10 2 sec, or even less, for m 1 / 2 and/or mo in the TeV range or for small 
gravitino mass mg. 

When the NLSPs decay, the high energy electromagnetic and hadronic particles are 
produced and may interact with background light nuclei and change the abundances. Using 
the output from low energy spectrum and interactions, we calculate the spectrum of decay 
products of NLSP and implement into BBN calculation to predict the primordial light 
element abundances. Then we require that the light element abundances after decay of 
NLSP are within the bounds of those inferred from the observational data. 

3.1 Collider and cosmological constraints 

We apply the same experimental bounds as in [14, 16]: (i) the lightest chargino mass 
m ± > 104 GeV; (ii) the lightest Higgs mass m h > 114.4 GeV; (iii) BR(B -> X s -f) = 
(3.55 ±0.68) x 10~ 4 and (iv) the stau mass bound m^ > 87 GeV. In this analysis we also 
use the top quark mass mj = 172.7 GeV [65]. However slight change of these experimental 
limit does not modify the final results, since the BBN constraints (especially from 6 Li/ 7 Li) 
are more severe. 

For the observational constraints on the primordial light element abundances, we use [6] 

1.2 x 10" 5 < D/H < 5.3 x 10~ 5 

Y p < 0.258 
8.5 x lO" 11 < 7 U/H 

SUe/D < 1.52 

^i/Li < 0.1 (0.66). 
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Figure 1: The plane (roi/2,mo) for tan/3 = 10, mg = too (left window) and tan/3 = 50, = 
0.2mo (right window) and for Ao = 0, /i > 0. The light brown regions labeled "LEP x + " and "LEP 
Higgs" are excluded by unsuccessful chargino and Higgs searches at LEP, respectively. In the right 
window the darker brown region labeled "& — > 37" is excluded assuming minimal flavor violation. 
The dark grey region below the dashed line is labeled "Tachyonic" because of some sfermion masses 
becoming tachyonic, and is also excluded. In the rest of the grey region (above the dashed line) the 
stau mass bound > 87 GeV is violated. In the region "No EWSB" the conditions of EWSB 
are not satisfied. The dotted diagonal line marks the boundary between the neutralino (x) and the 
stau (f) NLSP. The regions excluded by the various BBN constraints are denoted in violet. A solid 
magenta curve labeled "CMB" delineates the region (on the side of label) which is inconsistent with 
the CMB spectrum. In both windows, the dark green bands labeled "Tr = 10 7 GeV" and "10 8 " 
correspond to the total relic abundance of the gravitino (from the sum of thermal and non-thermal 
production), for a given (denoted) reheating temperature, lying in the favored range. In the light 
green regions (marked "NTP") the same is the case for the relic abundance from NTP process 
alone. The blue dashed line denotes the relaxed boundary of the BBN constrains when we use the 
conservative limit ^i/Td < 0.66. 

For "Li/T/i we will also use the conservative upper limit which is given in the bracket. It 
allows for the possibility of stellar ^i (and 7 Li) depletion. Since is more fragile than 
7 Li, post-BBN lithium processing may conceivably reduce the ^Li/TLi ratio. The reader is 
referred to Ref. [6] for a more detailed discussion of the adopted limits. In the first two 
figures below, the regions excluded by BBN constraints will be shaded violet and marked 
"BBN". We also show the BBN constraints with the conservative bound on ^i/li with 
blue dashed lines. 

As regards the total gravitino relic abundance Q^h 2 , we apply the 3 a range derived 
from WMAP 5 year data [66] 

0.091 < n d h 2 < 0.128, (3.1) 
which in the figures below will be marked as green bands and labeled "Q^h 2 " . 
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Figure 2: The same as fig. 1 but for a fixed gravitino mass, tan (3 — 10, mg = 10 GeV (left window) 
and tan/3 = 50, = 100 GeV (right window). 

As previously in [16], we also include the bound on the possible distortion in the nearly 
perfect black-body shape of the CMB spectrum [67] by the injection of energetic photons 
into the plasma. However we note that this constraint (delineated with magenta line with 
a label "CMB" over it) seems generally less important than that due to the BBN [68, 69]. 

3.2 Numerical results 

In figs. 1 and 2 we give the updated results with the constraints and the relic abundance 
of the gravitino in the same format as used in [16] to facilitate the comparison. Generally, 
the regions in white are left allowed after applying collider and BBN constraints. Green 
bands denote ranges of parameters where the total gravitino abundance reproduces the 
DM abundance (3.1), while in the light green bands (marked "NTP") only the NTP part 
of the gravitino abundance agrees with that range. 

In fig. 1 we present two cases with gravitino mass related to the soft scalar mass rriQ. 
In the left window we take tan/3 = 10 and mg = rriQ and in the right window tan/3 = 50 
and wig = 0.2m,Q. In fig. 2 we present two cases with a fixed mg. In the left window we 
fix tan (3 = 10 and mg = 10 GeV while in the right one tan (3 = 50 and mg = 100 GeV. 

We can see that, as before, the whole neutralino NLSP region (above the black dotted 
diagonal line) is again ruled out, as well as a part of the stau NLSP region corresponding 
to smaller (and therefore smaller mf and hence larger stau lifetimes). By comparing 
with the corresponding figures in [16], we can see that including the bound state effects 
causes the allowed region of stau NLSP consistent with DM abundance (green band) to be 
more constrained, and also strengthens an upper bound on Tr to 

T R < a few x 10 7 GeV. (3.2) 
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(For comparison, without the bound state effect, the bound was a few times 10 8 GeV [16].) 
This bound is consistent with the result of ref. [33, 34] where the bound state effect was 
mimicked by making a simple and approximate constraint on the lifetime of the stau, 
r < 5 x 10 3 sec. However, as we shall now show, in a number of cases the stau lifetime that 
can be consistent with avoiding strong bound-state catalytic effect can be almost an order 
of magnitude larger. Furthermore, by applying more conservative bounds on ^Li/^Li, the 
upper bound (3.2) can in some cases be violated by up to an order of magnitude. 

In order to examine the upper bound on Tr more closely, in fig. 3 we plot Tr vs. mg for 
all the points in our scans which satisfy the constraints from BBN and collider experiment 
and for tan/3 = 10 (left) and tan/3 = 50 (right). In these plots we use several choices of 
the gravitino mass: linear in the scalar mass (0.2mo, 0.4mo), linear in the gaugino mass 
(0.2m^ 2 , 0.4m 1 / 2 ) or constant (24 values between 0.1 GeV and 100 GeV). For each point, 
Tr is calculated using eqs. (2.2) and (2.6) for a given gravitino mass and for the mass spectra 
and the yield of NLSP at each point of parameter space, imposing n d h 2 = 0.11. Points 
marked by green dots survive all the constraints from experiments and cosmology while 
those marked by grey dots are allowed by experiments but disallowed by BBN constraints 
from light element abundances. Again, we can see that, for both choices of tan /3, we find 
the upper bound (3.2). This result does not change if one alters collider constraints since 
the bound is set by the 6 Li abundance. 

However, when we use the conservative bound on 6 Li/'Li < 0.66 then we find new 
isolated, albeit rather small allowed regions (marked by blue dots), where the reheating 
temperature can go up to 3 x 10 8 GeV. (They correspond to vertical blue dashed region 
around m 1( / 2 ~ 1 TeV and small ttiq in the right window in figs. 1 and 2). Clearly, those 
points evade the simple bound r < 5 x 10 3 sec (marked with a red dashed line). On the 
other hand, this region corresponds to a rather low stau mass around 100 GeV (but still 
above the LEP limit). This can be seen in fig. 4 where we plot stau lifetime for the data 
points shown in fig. 3. We can see that the (green) points allowed by 6 Li/ 7 Li < 0.1 are 
located below the stau lifetime of 5 x 10 3 sec for tan (3 = 10 and 3 x 10 4 sec for tan (3 = 50. It 
is also clear that, at smaller tan (3 (left panel) and larger rrif the bound r < 5 x 10 3 sec can 
actually be too weak. On the other hand, assuming the conservative bound 6 Li/ 7 Li < 0.66 
(blue points), at stau mass around 100 GeV, for tan (3 = 50 the lifetime can exceed 10 4 sec 
even by a few orders of magnitude. Such low stau mass will be accessible to the LHC, thus 
allowing one to scrutinize such cases fairly easily. 

The reason why the points corresponding to stau NLSP with such a long lifetime are 
allowed is that, as the stau mass decreases, its yield Y also decreases and can drop below 
10~ 14 . In this case the BBN constraints from the bound state effects can be avoided even 
with stau lifetime as long as 3 x 10 4 sec (green points). 

To see this, we present the left window of fig. 5. Assuming a standard (more conser- 
vative) limit on ^Li/TLi we find green (blue) points corresponding to the lifetime of up to 
3x 10 4 sec (2x 10 7 sec) and the stau yield of less than 10~ 14 . For lifetimes longer than 10 7 sec 
the constraint due to violating the upper limit on 3 He/D due to 4 He photo-disintegration 
becomes severe [31]. On the other hand, for much longer lifetimes, between 10 10-11 sec, we 
again find some isolated allowed regions of (blue) points. They become allowed because 
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Figure 3: The Tr vs. plane for tan/3 = 10 (left) and tan/3 = 50 (right). Data points 

which survived all BBN and experimental constraints are shown with green dots, while grey dots 
are disallowed by the BBN while allowed by collider constraints. Blue points are added to the 
green when we use a more conservative limit 6 Li/ 7 Li < 0.66. Tr is determined so that the total 
(thermal and non-thermal) production of gravitinos satisfies the WMAP value; here we used il^h 2 = 
flj^h 2 + £l~ TP /i 2 = 0.11. We also show the lifetime contour of stau, r = 5 x 10 3 sec, which 
was obtained using the specific relation between stau mass and gaugino masses as well as the 
approximate relic density of the f in the CMSSM. 

the mass difference between the gravitino LSP and the stau NLSP is very small there and 
the electromagnetic energy released is not enough to be a potential problem for BBN. 

Finally, in the right window of fig. 5 we show Tr vs. stau lifetime for all the cases we 
have considered, and which summarizes many of our points made above. In particular, with 
the ranges of primordial light element abundances adopted in this analysis (green points), 
the upper limit on Tr scales with the stau lifetime t but does not exceed the limit (3.2). 
Allowing the more conservative limit on "Li/ "Li pushes it up to nearly 10 8 GeV, with some 
isolated points tolerating Tr as high as 3 x 10 8 GeV, as stated above. 

4. Summary 

We have re-analyzed the gravitino as dark matter in the Universe in the framework of the 
CMSSM, taking into account of a number of recent improvements in the calculation. The 
improvements concern the thermal production of gravitinos, including the decay allowed 
by thermal masses as well as scatterings, the computation of the hadronic spectrum from 
NLSP decay including the correct implementation of left-right stau mixing, as well as an 
updated BBN calculation fully treating the effects of bound state effects between charged 
massive particles and nuclei. We found that the over-abundance of 6 Li from bound state 
effects puts more severe constraints on allowed ranges of CMSSM parameters than before 
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Figure 4: The same as fig. 3 but in the plane of lifetime and mass of stau for tan/3 = 10 (left) and 
tan/3 = 50 (right). 
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Figure 5: The same as fig. 3 but in the plane of the stau yield Y = n/s (left window) and of Tr 
(right window) vs. stau lifetime for both tan/3 = 10 and tan/3 = 50. 

and puts an upper bound on Tr at a few times 10 7 GeV and the stau lifetime less than 
around 3 x 10 4 sec. 

We also analyzed the impact of applying a more conservative upper limit on 6 Li/ 7 Li. In 
this case a new region at small stau mass opens up because of the reduced relic abundance 
of the stau, even though the lifetime of the stau is much longer than 10 4 sec, and can reach 
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up to 5 x 10 sec, in some cases even between 10 and 10 sec. In these latter cases, which 
are easily testable at the LHC, Tr can exceed 10 8 GeV. 
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